Electron-Phonon Interactions for Optical Phonon Modes in Few-Layer Graphene 
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We present a first-principles study of the electron-phonon (e-ph) interactions and their contribu- 
tions to the linewidths for the optical phonon modes at F and K in one to three-layer graphene. It 
is found that due to the interlayer coupling and the stacking geometry, the high-frequency optical 
phonon modes in few-layer graphene couple with different valence and conduction bands, giving 
rise to different e-ph interaction strengths for these modes. Some of the multilayer optical modes 
derived from the F-E2 g mode of monolayer graphene exhibit slightly higher frequencies and much 
reduced linewidths. In addition, the linewidths of K-A[ related modes in multilayers depend on the 
stacking pattern and decrease with increasing layer numbers. 
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The possibilities of developing carbon-based nanos- 
tructures for electronics applications have stimulated re- 
cent interest in graphene and its derivatives^ 2 -^. One of 
the focus areas is to understand the scattering processes 
of electrons. In carbon nanotubes and graphite, the high 
currents or optical excitations have been shown to induce 
a significant overpopulation of the optical phonon modes 
of E 2g at T (T-E 2g ) and A[ at K (K-A[)&&&. Since these 
phonon modes exhibit a strong electron-phonon (e-ph) 
interactions, overpopulation of them leads to a dramatic 
reduction of the ballistic conductance of carbon nan- 
otubes at high bias potentials^, and consequently inter- 
connect performance deteriorates^. Understanding the 
phonon decays 8 from a microscopic point of view, in par- 
ticular, based on the e-ph interaction, is thus a key step to 
improve the transport properties of these materials and 
to control device performance. Furthermore, the e-ph in- 
teraction also plays a significant role in many phenomena 
such as the quasiparticle dynamics and anomalies in pho- 
toemission spectra 9 - ' 10 ? 11 ! 12 ? 13 , Raman scattering^, and 
superconductivity^. 

Experimentally, the linewidths of the zone-center 
phonon modes obtained from Raman or infared (IR) 
measurements contain significant contributions from the 
e-ph interaction 16 . In graphene, the phonon linewidth 
of the T-E 2g mode is estimated to be about 13 cm -1 
based on the Raman spectra^. In graphite, the graphene 
E 2g mode splits into two branches: the Raman-active 
E 2g and IR-active E\ u modes. The linewidth of the 
Raman-active mode (11.5 cm" 1 —) is almost the same as 
that of graphene, while IR measurements show that the 
linewidth of the E\ u mode is surprisingly much smaller 19 . 

Few-layer graphene (FLG) presents an interesting sys- 
tem because of the possibility to tune its electronic 
properties^ .. In the epitaxially grown graphene, FLG 
is often produced as a main product* 2 -. Depending on 
the layer number and stacking geometry^, the linear 
band dispersions in monolayer graphene evolve into sev- 
eral bands due to the interlayer coupling in FL G 21 ' 22 . 
Similarly, the E 2g mode at T and A\ mode at K in sin- 
gle layer also split into several branches^. It is expected 
that the e-ph interaction will be significantly modified as 
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FIG. 1: (a) Supercell (shaded area) used in the frozen-phonon 
calculation for the K phonon modes in few-layer graphene. 
(b) First Brillouin zone of the primitive unit cell of graphene. 
The square area for the dense k-grid sampling is indicated. 



the number of layers and stacking geometry are varied. 

In this work, we performed first-principles calculations 
of the e-ph interactions in one-, two- (AB stacking), and 
three-layer (ABA and ABC stackings) graphene. The T- 
E 2g and K-A[ modes in graphene, which are respectively 
coupled to intra- and inter-valley electronic scatterings, 
have much stronger e-ph interactions than other modes 2 ^. 
Here we focus on these two modes and their derivatives 
in FLGs in order to investigate the stacking effect. We 
found that while the weak interlayer interaction gives rise 
to only small splittings of these phonon modes, the result- 
ing e-ph interactions for some of them are considerably 
suppressed because of symmetry constraints. 

The e-ph matrix element Soc+qW kj 1S defined as 
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where SV sc f = K, c /(iiq) — V sc f(0) is the variation of the 
self-consistent potential field due to the perturbation of 
a phonon with wave vector q and branch index v. |k, j) 
is the electronic Bloch state. We further define the e-ph 
coupling strength between phonon mode qy and Bloch 
states |k + q,j') and \k,j) as ^-(k) = \g" k+cl)f tWj \ 2 ■ 

The electronic states and the self-consistent field 
were computed using the first-principles codes Quan- 
tum ESPRESSO 2 ^! 2 ^ and the perturbation of a phonon 
mode was handled with the frozen-phonon approach. 
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FIG. 2: (Color online) (a)-(d) Optical modes at F, (e)-(h) band dispersions near the Fermi level, and (i)-(l) the square of 
the e-ph interaction strength |g| 2 of the optical modes at V for monolayer, AB bilayer, ABA and ABC trilayer graphene, 
respectively. Symmetry-allowed transitions from valence bands to conduction bands with nonzero strength are shown. 



Fig. [Ha) shows the supercell used in the frozen-phonon 
calculation for the K phonons. The electronic struc- 
ture was calculated with the local density approxima- 
tion (LDA) within density-functional theory and the 
core- valence interaction was modeled by norm-conserving 
pseudopotentials 2 ^. Wave functions of the valence elec- 
trons were expanded in plane waves with a kinetic energy 
cutoff of 70 Ry. The phonon frequencies and associated 
eigenvectors were computed using the density-functional 
perturbation theory (DFPT) 25 , details of which have 
been presented in our previous workS. A vacuum re- 
gion of 10 A was introduced in our supercell to eliminate 
the artificial interaction between neighboring supercells 
along the z direction. The relaxed C-C bond length is 
1.42 A and the interlayer distance is 3.32 A for all FLGs 
considered in this paper. Variations of the potential fields 
SV sc f were calculated through self-consistent calculations 
to find the potential field for both perturbed and unper- 
turbed systems. The following calculations were carried 
out on a dense 100x100 fc-grid within a small square 
area enclosing point K in reciprocal space, as indicated 
in Fig. [IJb) in order to obtain the electronic wave func- 
tions at k and k + q near the Fermi level. This dense 
/c-sampling was found necessary for a quantitative de- 



scription of the scattering process near the Fermi level in 
one- and few-layer graphene. Finally the e-ph interac- 
tion matrix elements were computed using Eq. (1). 

In order to check the accuracy of our calcula- 
tions, we first calculate the e-ph matrix elements 
over the Fermi surface and compare them with 
previously published results for monolayer graphene. 
Due to the electronic degeneracy at K, the aver- 
aged e-ph matrix elements for all possible pairs are 
(Sr>F=£i. 3 ; I <7(K)i,Ki 174=0.0401 eV 2 for the T-E 2g 
mode, and (fl&>F=Eij ■ l<7(2K)i,K,-| 2 /4=0.0986 eV 2 for 
the K-A[ mode, respectively. These results are in excel- 
lent agreement with those in previous DFPT calculations 
(0.0405 and 0.0994 eV 2 , respectively)^. 

The phonon linewidth 7 due to the e-ph coupling is 
defined as 2 ^: 

x5[e ki -e (k+qW , (2) 

with /iy being the Fermi-Dirac occupation function for 
Bloch state |k,j). We used a broadening parameter of 
0.01 eV for the <5-function in Eq. ©, and k B T=2.5 meV 
for the Fermi-Dirac distribution. 
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FIG. 3: (Color online) Schematic diagrams of electronic or- 
bitals and optical phonon modes at T that give rise to non-zero 
e-ph matrix elements in bilayer graphene. The tight-binding 
orbitals of each band close to K are shown. Red (blue) is 
positive (negative) amplitude. 



Table Q] lists the calculated phonon linewidths due to 
the e-ph coupling. For the monolayer, the linewidth 7 
of the degenerate T-i?2 ff modes is 11.2 cm -1 , in rea- 
sonable agreement with experimental observation of 15.0 



cm — and previously published result of 11.5 cm' 
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A value of 20.4 cm" 1 is obtained for the highest op- 
tical K-A[ mode. In FLGs, the phonon modes split 
into branches with different symmetries. Consequently, 
their linewidths can be quite different from one another. 
In the AB bilayer, for example, the original highest T 
phonon splits into the Raman active E g and IR-active 
E u modes. While the linewidth of the E g mode (8.6 
cm -1 ) is only slightly smaller than that of the monolayer 
(11.2 cm -1 ), the linewidth of the E u mode is two orders 
of magnitude smaller. The calculated linewidth of the E g 

I 



mode is consistent with recent experimental observation 
of 13.5 cm -1 for the V optical phonon mode in bilayer 
graphene^. Similar results for the monolayer and bi- 
layer graphene were also reported by Park et al£^. For 
the ABA and ABC trilayer, the results in Table Q] show 
different behavior for different modes. In particular, for 
the ABC trilayer two out of three optical modes at V 
show reduced electron-phonon interaction. Most impor- 
tantly, the linewidths of all K phonons in Table Q] drop 
significantly as the layer number increases from 1 to 3. 
This indicates that the interlayer interaction can effec- 
tively suppress vallley-spin decoherence via e-ph scatter- 
ings, making the FLGs more attractive as valleytronic 
materials22i22. 

For FLGs there are a few valence and conduction 
7T bands near the Fermi energy level. The symmetry 
allowed interband transitions by the absorption of a 
phonon are indicated in the middle column of Fig. [Jl To 
examine the contributions of different electronic states to 
the e-ph coupling strengths in FLGs, we present in the 
right column of Fig. [J] the absolute value squared of the 
e-ph coupling matrix elements for optical phonons at T 
as a function of electronic crystal momentum k along the 
symmetry line T-K for all symmetry allowed transitions. 
Some of these matrix elements decrease monotonically 
with increasing k x , as in the case of monolayer; some of 
them show a minimum at K. These features are closely 
related to the symmetry of the electronic states near K in 
the presence of interlayer interactions, and can be further 
quantitatively understood using a tight-binding model. 
Following the procedures in Refs. [34] and [HI, we obtain 
the e-ph matrix element for the in-plane optical phonon 
mode qv in L-layer graphene as: 
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with C/t(k)_ = [ul a:j (k),ul p j{k),u^(k),u* 0J (k),...} 
being the tight-binding amplitudes of band j for each 
site in the unit cell. e} Q and tip are the vibrational eigen- 
vectors for the two atoms in layer /. f(k) = Yli=i ^e lk Ri , 
where <5,'s (z=l-3) are the unit vectors connecting atom 
a in layer I to its three nearest- neighbors (NNs) , and R; 
are the lattice vectors of the unit cells in which the three 
NNs are located. The constant g = JQ 1 / 2 / '(oj^s/M) de- 
pends on the mode frequency co", the e-ph interaction 
parameter J, the area of the unit cell f2, and the carbon 



atomic mass M. The coupling matrix < & l/ (k, q) has the 
form: 



*"(k,q) 



with 
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TABLE I: Calculated phonon linewidth 7 (in cm 1 ) for the high-frequency optical phonon modes at V and K in monolayer, 
bilayer, and trilayer graphene. The mode symmetries S and the frequencies lo (in cm -1 ) are also listed for completeness. 







Monolayer 






AB 






ABA 






ABC 




s 




7 


S 


Id 


7 


S 




7 


S 


Id 


7 


r 


E 2g 


1586 


11.2 


E 3 


1587 


8.6 


K 


1586 


9.7 


Eg.a 


1586 


7.2 










E u 


1592 


0.1 


E" 


1588 


11.0 


E u 


1589 


0.0 
















K 


1593 


2.8 


E g,b 


1594 


0.3 


K 




1306 


20.4 


E 


1318 


9.0 


E[,E'l 


1316 


8.4 


E 


1318 


2.8 
















E' 2 


1324 


3.6 




1325 


2.2 



t(k) 



(q) 



r 



^(q)-*( k - 




q)-^(q) 



(6) 



(a)E 



2g 



(b) E „ 





FIG. 4: (Color online) Electronic momentum-resolved contri- 
butions of all scattering processes to the phonon linewidth for 
the doubly degenerate (a) E 2g mode in monolayer graphene 
and (b) E g mode in bilayer graphene. Note that only - 
O.OOg^fc^fcy^O.OOg for monolayer and -0.021<fc :E ,fc a <0.021 
for bilayer are shown. k x and k y are in 2n/ao- 



It is clear that the e-ph interaction strength depends 
on the orbital characteristics of the initial and final elec- 
tronic states as well as the displacement pattern of the 
phonon mode. Note that is an off-diagonal matrix, 
coupling the electronic component of one sublattice to 
that of the other sublattice in each layer through phonon 
displacements in the same layer. The coupling effects in 
different layers will be summed up in either a constructive 
or a destructive way, depending on the relative phase of 
the phonon displacements in different layers and the rel- 
ative phase of the tight-binding amplitudes in the initial 
and final electronic states involved. 

As an example, Figure [3] schematically shows the 
symmetry-allowed transitions through the E g and E u 
modes at T between valence and conduction bands in 
bilayer graphene. The four relevant energy bands shown 
in Fig. [5{f) are indexed as 2, 1, 1', and 2' from lower 
valence bands to higher conduction bands. For the E g 
mode, the in-plane vibrations of a and (3 atoms in neigh- 
boring planes are in phase. Taking into account of the 
symmetry of electronic wave function in each layer, al- 
lowed transitions turn out to be 1 — > 1' and 2 — > 2'. As k 
approaches K, the four bands become parabolic and the 



wave- function amplitudes tend to be localized on one sub- 
lattice of each layer. Since the transition matrix elements 
for 1— ► 1' and 2— > 2' contain product of the electronic 
wave functions on both sublattices, small wave function 
amplitudes give rise to vanishingly small matrix elements 
near K, as shown in Fig. [Uj). On the other hand, the 
E u mode involves in-plane vibrations of a and (3 atoms 
in neighboring planes that are out of phase. The sym- 
metry allowed transitions are 2 — > 1' and 1 — > 2'. Even 
though the wave function amplitude at one sublattice 
may become small when k approaches K, the products 
with wave function amplitudes at the other sublattice are 
still finite. Therefore, the matrix elements exhibit similar 
linear behavior as in monolayer graphene. 

In collecting the contributions from different k states 
to calculate the phonon linewidth in Eq. ([2]) , energy con- 
servation is controlled by the delta function. For the E g 
mode, this condition can be satisfied in the transition 
1 — > V but not in the transition 2 — > 2'. Neverthe- 
less, the phonon linewidth listed in Table I is about 3/4 
of the value for monolayer graphene. Figure 2] presents 
the k-resolved linewidth contributions for the doubly de- 
generate E2 g mode in monolayer graphene and the E g 
mode in bilayer graphene. Clearly, the contributions are 
almost isotropic in k space for the Ei g mode in mono- 
layer graphene. The shell shape is a result of energy 
conservation for the scattering events. Due to the trig- 
onal symmetry of the constant-energy surface in bilayer 
graphene, the k-dependent contributions show a three- 
fold symmetry for the E g mode. In contrast, the symme- 
try allowed transitions of 2 — > 1' and 1 — > 2' for the E u 
mode cannot satisfy the required energy conservation be- 
cause the energy separation of the initial and final states 
is too large. This leads to a vanishingly small probability 
for the E u mode to decay through the e-ph interaction. 
These results are a reminiscence of the properties of T- 
E 2g and T-Ei u phonons in bulk graphite regarding their 
e-ph interactions^. 
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FIG. 5: (Color online) Electronic momentum-resolved contri- 
butions of all scattering processes to the phonon linewidth for 
the (a) E' a , (b) E", and (c) E' b phonon modes at P in ABA 
trilayer graphene. The k area shown in the figure is indicated 
in Fig. [ljb). The origin corresponds to K. 

For the ABA trilayer, the four energy bands crossing 
or touching the Fermi energy are again indexed as 2, 1, 1' 
and 2' in Fig. [2] When a T-E' a phonon is absorbed, both 
symmetry and energy allowed transitions are 2— >2' and 
1— >1', which manifest themselves as the inner and outer 
rings respectively in Fig. EJa) . When a Y-E" phonon is 
absorbed, the symmetry and energy allowed transitions 
are 2^ V and 1— » 2'. Incidentally, the energy differences 
between bands 2 and 1' and between bands 1 and 2' are 
equal for any k. Therefore, the rings associated with 
these two different transitions overlap in Fig. [5]Jb) . The 
electronic states involved span an energy range of Hlj ~ 
0.2 eV about the Fermi level. A Y-E' b phonon obeys the 
same selection rule as a Y-E' a phonon since both bear the 
same group representation. The e-ph coupling strength 
of the former is, however, much weaker than that of the 



latter. 

For the ABC trilayer, the relevant energy bands are in- 
dexed as 1 and 1' in Fig. [21 The three Y phonons derived 
from the Ei g mode in monolayer graphene are respec- 
tively denoted by T-E g>a , Y-E u , and Y-E g ^ in the order 
of increasing frequencies. Interestingly, only the Y-E ga 
phonons show an appreciable e-ph coupling correspond- 
ing to the transition of 1 — * 1'. 

In summary, we have studied the phonon linewidths 
of the high-frequency optical phonon modes in few-layer 
graphene. We found that there is a strong suppression 
of the e-ph interaction for these modes resulting from 
stacking patterns. The split optical phonon modes in 
few-layer graphene are shown to only couple with the 
electronic bands of specific orbital symmetry and exhibit 
various e-ph interaction strengths. These features are 
well illustrated using a tight-binding model. 
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